Rheology of Soft Glassy Materials 
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We attribute similarities in the rheology of many soft materials (foams, emulsions, slurries, etc.) 
to the shared features of structural disorder and metastability. A generic model for the mesoscopic 
dynamics of "soft glassy matter" is introduced, with interactions represented by a mean-field noise 
temperature x. We find power law fluid behavior either with (x < 1) or without (1 < x < 2) a yield 
stress. For 1 < x < 2, both storage and loss modulus vary with frequency as uj x ~ x , becoming flat 
near a glass transition (x — 1). Values of x ~ 1 may result from marginal dynamics as seen in some 
spin glass models. 
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Many soft materials, such as foams, emulsions, pastes 
and slurries, have intriguing rheological properties. Ex- 
perimentally, there is a well-developed phenomenology 
for such systems: their nonlinear flow behavior is often 
fit to the form a = A + Bj n where a is shear stress 
and 7 strain rate. This is the Herschel-Bulkeley equa- 
tion f||; or ( for -4 = 0) the "power-law fluid" |]-|. 
For the same materials, linear or quasi-linear viscoelastic 
mesurements often reveal storage and loss moduli G'(ui), 
G"(u>) in nearly constant ratio (G" /G' is usually about 
0.1) with a frequency dependence that is either a weak 
power law (clay slurries, paints, microgels) or negligible 
(tomato paste, dense emulsions, dense multilayer vesicles, 
colloidal glasses) [Q-[L(i. This behavior persists down to 
the lowest accessible frequencies (about 10~ 3 -1 Hz de- 
pending on the system), in apparent contradiction to 
linear response theory O], which requires that G"{ui) 
should be an odd function of to. 

That similar anomalous rheology should be seen in 
such a wide range of soft materials suggests a common 
cause. Indeed, the frequency dependence indicated above 
points strongly to the generic presence of slow "glassy" 
dynamics persisting to arbitrarily small frequencies. This 
feature is found in several other contexts jl2H14| , such as 
elastic manifold dynamics in random media p^ , p^ ] . The 
latter is suggestive of rheology: charge density waves, 
vortices, contact lines, etc. can "flow" in response to an 
imposed "stress" . In this Letter we argue that glassy dy- 
namics is a natural consequence of two properties shared 
by all the soft materials mentioned above: structural 
disorder and metastability. In such materials, thermal 
motion alone is not enough to achieve complete struc- 
tural relaxation. The system has to cross energy bar- 
riers (for example those associated with rearrangement 
of droplets in an emulsion) that are very large com- 
pared to typical thermal energies. Therefore the sys- 
tem adopts a disordered, metastable configuration even 
when (as in a monodisperse emulsion or foam) the state 
of least free energy would be ordered ^] . While the im- 
portance of disorder has been noted before for specific 
systems |^, |TT| , |l8| -^lf , we feel that its unifying role in rhe- 
ological modelling has not been appreciated. 



To test these ideas, we construct a minimal "generic 
model" for soft glassy matter. For simplicity, we ignore 
tensorial aspects, restricting our analysis to simple shear 
strains. Consider first the behavior of a foam or dense 
emulsion under shear. We focus on a mesoscopic region, 
large enough for a local strain variable I to be defined, 
but small enough for this to be approximately uniform 
within the region, whose size we choose as the unit of 
length. As the system is sheared, droplets in this region 
will first deform elastically from a local equilibrium con- 
figuration, giving rise to a stored elastic energy (due to 
surface tension, in this example ps|). This continues up 
to a yield point, characterized by a strain l y , whereupon 
the droplets rearrange to new positions in which they 
are less deformed, thus relaxing stress. The mesoscopic 
strain I measured from the nearest equilibrium position 
(i.e., the one which can be reached by purely elastic de- 
formation) therefore executes a saw-tooth motion as the 
macroscopic strain 7 is increased p2[ . Neglecting nonlin- 
earities before yielding, the local shear stress is given by 
kl, with k an elastic constant; the yield point defines a 
maximal elastic energy E = ^kl y 2 . A similar description 
obviously extends to many others of the soft materials 
discussed above. 

We now ascribe to each mesoscopic region not only 
its own strain variable I, but also its own maximal yield 
elastic energy, E > 0. We model the effects of structural 
disorder by assuming a distribution of such yield ener- 
gies E, rather than a single value common to all regions. 
The state of a macroscopic sample is then characterized 
by a probability distribution P(l,E;t). We propose the 
following dynamics for the time evolution of P: 

§t p = ^§i p " ^~ (E ~ hkl2),x p + r(t) p(E)S(i) (i) 

The first term on the r.h.s. arises from the elastic de- 
formation of the regions. This embodies a mean field 
assumption, that between successive local yield events, 
changes in local strain follow those of the macroscopic 
deformation: I = 7. Note, however, that due to stochas- 
tic yielding events the stress kl is spatially inhomoge- 
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neous (as is the local strain I). The macroscopic stress is 
defined as an average over regions 



a(t) = k (I) = k / I P(l, E; t) dl dE 



(2) 



The second term on the r.h.s. of (Q) describes the yielding 
of our mesoscopic regions. We have written the yielding 
rate as the product of an "attempt frequency" Tq, and 
an exponential probability for activation over an energy 
barrier E — ^kl 2 (the excess of the yield energy over that 
stored elastically) . However, the resemblance to thermal 
activation is formal: we expect these "activated" yield 
processes to arise primarily by coupling to structural re- 
arrangements elsewhere in the system. In a mean-field 
spirit, all such interactions between regions are subsumed 
into an effective "noise temperature" , x. We first regard 
To, x as arbitrary constants, but later discuss their mean- 
ing and their possible dependences on other quantities. 

Finally, the third term on the r.h.s. of (|l|) describes the 
relaxation of regions to new local equilibrium positions 
after yielding, which we treat as effectively instantaneous. 
The first factor in this term is simply the total yielding 
rate T(t) — T (exp[— [E — ^kl 2 )/x]) p . The remaining 
two factors incorporate further mean-field assumptions, 
as follows. First, the yield energy E for distortions about 
any equilibrium configuration is uncorrelated with the 
previous one for this region; it is drawn randomly from 
the prior distribution ("density of states") p(E) which 
we assume to be time-independent. Second, immediately 
after yielding, a region always finds itself in a completely 
unstressed state of local equilibrium with I — (hence 
the Dirac delta function, 6(1)). This latter simplification 
is not essential, as shown elsewhere [ p3| . 

In the absence of flow (j(t) = 0), the model ([!]) 
describes activated hopping between "traps" of depth 
E' = E — \kl 2 with density p(E'). This corresponds 
to Bouchaud's model for glassy dynamics p^-p^] , whose 
predictions we briefly recall. For high (noise) tempera- 
tures x the system evolves towards the Boltzmann dis- 
tribution P eq (E') ~ p(E') exp(E'/x). As x is lowered, 
this distribution may cease to be normalizable, leading 



to a glass transition at x 



lim f 



y (d/dE) \np(E). 



For x < x g , no equilibrium state exists, and the sys- 
tem shows "weak ergodicity breaking" and various aging 
phenomena. A finite value of x g implies an exponen- 
tial tail in the density of states, p ~ exp(— E/x g ), which 
corresponds to a Gaussian distribution of yield strains 

ly = (2E/k) 1 / 2 . 

A major attraction of the model defined by (|l|) and (||) 
is that an exact constitutive equation, relating the stress 
a(t) to the strain-rate history [j(t' < t)], can be ob- 
tained Since this is quite complicated, we restrict 
ourselves here to two standard rheological tests, for which 
the full form is not required. We use non-dimensional 
units for time and energy by setting Tq = x g = 1; we 
also rescale our strain variables (l,j) so that k = 1. In 
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FIG. 1. Linear moduli G' (solid line) and G" (dashed) vs 
frequency uj at various noise temperatures. 

these units, p(E) = exp[-E(l+f(E))] with f(E) -> for 
E — > co. Up to sub-power-law factors such as logarithms, 
all power laws reported below are valid for any f(E); nu- 
merical examples use / = 0. Analytical and numerical 
support for our results will be detailed elsewhere |{23| . 

Consider first the complex dynamic shear modulus 
G*(oj) = G' + iG", which describes the stress response to 
small shear strain perturbations around the equilibrium 
state. As such, it is well defined (i.e., time-independent) 
only above the glass transition, x > 1. Expanding (^) 
to first order in the amplitude 7 of an oscillatory strain 
j(t) = -fcoaut, we find G*(u>) = (iuir/ (iuir + l)) eq . This 
corresponds to a distribution of Maxwell modes whose 
spectrum of relaxation times r = exp(E/x) is given by 
the equilibrium distribution P cq (E) ~ exp(E / x) p(E) . 
The relaxation time spectrum thus exhibits power law 
behavior for large r: P(r) ~ t~ x . This leads to power 
laws for G* in the low frequency range (Fig. [I]): 



G" 
G' 



for 2 < x, 
for 3 < x. 



for 1 < x < 2 
for 1 < x < 3 



(3) 



For x > 3 the system is Maxwell- like at low frequencies, 
whereas for 2 < x < 3 there is an anomalous power law 
in the elastic modulus. Most interesting is the regime 
1 < x < 2, where G' and G" have constant ratio; both 
vary as u> x ~ 1 . Behavior like this is observed in a number 
of soft materials f7|JTc|l . Moreover, the frequency expo- 
nent approaches zero as a: — > 1, resulting in essentially 
constant values of G" and G', as reported in dense emul- 
sions, foams, and onion phases [^|-^|. Note, however, that 
the ratio G" /G' ~ x — 1 becomes small as the glass tran- 
sition is approached. This increasing dominance of the 
elastic response G' prefigures the onset of a yield stress 
for x < 1 (discussed below) If a high energy cutoff 
B max is imposed on p(E) (giving an upper limit on lo- 
cal yield strains) , the above results remain valid down to 
^min = exp(— E max /x). Well-defined equilibrium values 
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FIG. 2. Shear stress a vs shear rate 7, for a; = 0.25, 0.5, 
. . ., 2.5 (top to bottom on left); x — 1, 2 are shown in bold. 
Inset: small 7 behavior, with yield stresses for x < 1 shown 
by arrows. 

of the linear moduli then exist also for x < 1; one still 
finds G" ~ w 31-1 for w m i n < w < 1. For x just below 
x g = 1, a log-log plot of G"(uj) therefore exhibits a small 
negative slope (whereas G" is constant). This may again 
be compatible with recent experimental data J^-|l0|. 

We now turn to the case of steady shear flow, 7 = 
const, for which the steady state distribution P ss {l,E) 
can be obtained analytically. After integrating over E, 
one finds P ss (l) ~ Q(l)g(z(l)) with 

z(l) = ~ f Q el2/2xd l 9( z ) = J P( E ) exp(-ze~ B/x ) dE 

In the large z limit, g(z) ~ z~ x . Figure ^ shows that 
for large shear rates, 7 > 1, a increases very slowly for 
all x (a ~ (xlnj) 1 / 2 ). More interesting is the small 7 
behavior, where we find three regimes: (i) For x > 2, the 
system is Newtonian, a = 777. The viscosity is simply the 
average relaxation time 77 = (exp(E / x)) eq = (t) taken 
over the equilibrium distribution of energies, P eq (E) ~ 
exp(E/x)p(E). Hence 77 ~ (exp(2E/x)) , which diverges 
at x — 2. (ii) For 1 < x < 2 one finds power law fluid 
behavior, a ~ -y^ -1 . (iii) For x < 1, the system shows a 
yield stress: 17(7 — > 0) = cr y > 0. (This has a linear on- 
set near the glass transition, <j y ~ 1 — x.) Beyond yield, 
the stress again increases as a power law of shear rate, 
cr — <7 y ~ 7 1 ~ a: (for 7 <C 1). The behavior of our model in 
regimes (ii) and (iii) therefore matches respectively the 
power-law fluid [l)-|j| and Herschel-Bulkeley [Q|2) sce- 
narios as used to fit the nonlinear rhcology of pastes, 
emulsions, slurries, etc. 

We now speculate on the origin and magnitude of the 
"attempt frequency" Tq and the "noise temperature" 
x. First note that the parameter Tq is the only source 
of a characteristic timescale (chosen as the time unit 
above). We have approximated it by a constant value: 



To (7) = To(0). One possibility is that the intrinsic rate 
constant Tq arises from true thermal processes. If so it 
can be estimated as T\ oc kBTP oq (0)Q with ri oc a local 
diffusive attempt rate (for 1 /xm emulsions this might 
be 0.01s); fcBTP O q(0) is the (small) fraction of regions 
in which true thermal activation can surmount the yield 
barrier. The factor Q denotes the number of neighboring 
regions perturbed as a result of one such thermal event. 
A more detailed analysis (involving an extension to our 
model ^3|) then shows that kBTP cq (Q)Q must be large 
enough (at least of order unity) to avoid depletion of the 
low energy part (E < k^T) part of the barrier distri- 
bution. This mechanism may arise in systems (such as 
foams) in which one local rearrangement can trigger a 
long sequence of others pO|pl| . If so, the resulting in- 
trinsic rate Tq ~ ri oc provides a plausible rheological 
timescale. (If Q is too small, Tq will instead be of order 
F\ oc e~ E ^ kBT , which for typical barrier energies E = (E) 
is unfeasibly slow.) 

We emphasize, however, that Tq may be strongly sys- 
tem dependent, and any specific interpretation of it re- 
mains speculative. Nonetheless we may view the acti- 
vation factor in Eq.(l) as the probability that a pertur- 
bative "kick" to a given mesoscopic region (from events 
elsewhere) causes it to yield. We believe this activation 
factor should be primarily geometric in origin and hence 
depend on the disorder, but not on any intrinsic energy 
scale. Accordingly (in our units) x values generically of 
order unity can be expected. We argue next that x values 
close to unity may be normal. 

Consider first a steady shear experiment. For soft 
metastable materials, the rheological properties of a sam- 
ple freshly loaded into a rheometer are usually not repro- 
ducible; they become so only after a period of shearing 
to eliminate memory of the loading procedure. In the 
process of loading one expects a large degree of disor- 
der to be introduced; the initial dynamics under flow 
should therefore involve a high noise temperature x S> 1 . 
As the sample approaches the steady state, the flow will 
(in many cases) tend to eliminate much of this macro- 
scopic disorder pj| so that x will decrease. But, as this 
occurs, the noise-activated processes will slow down; as 
x — > 1, they become negligible. Assuming that, in their 
absence, the disorder cannot be reduced further, x is then 
"pinned" at a steady-state value at or close to the glass 
transition. This scenario, although extremely specula- 
tive, is strongly reminiscent of the "marginal dynamics" 
seen in some mean-field spin glass models [^6| . 

There remains several ambiguities within this picture, 
for example whether the steady state value of x should 
depend on 7; if it does so strongly, our results for steady 
flow curves will of course be changed. If a steady flow 
is stopped and a linear viscoelastic spectrum measured, 
the behavior observed should presumably pertain to the x 
characterizing the preceding steady flow (assuming that x 
reflects structure only). But unless the strain amplitude 
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is extremely small the x-value obtained in steady state 
could be affected by the oscillatory flow itself |27| . 

Also uncertain is to what extent a steady energy input 
is needed to sustain the nonlinear dynamics. Although 
not represented in the model, a small finite strain rate 
amplitude might be needed to balance the gradual dis- 
sipation of energy in yield events. In its absence, one 
might expect the sample to show aging (i.e., P(l,E;t) 
nonstationary in time). Within the model, aging in fact 
only occurs for x < 1 (the regime for which we pre- 
dict a yield stress). Conversely we saw above that, even 
in this regime, for finite 7 a well-defined steady state dis- 
tribution is recovered: flow interrupts aging ]i"3| . This 
can be understood by considering the distribution of en- 
ergies. Without flow, one obtains a Boltzmann distri- 
bution P(E) ~ p(E)e E/x up to (for x < 1) a cutoff 
which shifts to higher and higher energies as the system 
ages |l^] . This cutoff, and hence the most long-lived traps 
visited (which have a lifetime comparable to the age of 
the system) , dominate the aging behavior Q. The pres- 
ence of flow leads to a steady state value of this cutoff 
of E ~ x\n( A /~ 1 x 1 / 2 ), while for higher energies one has 
P SS (E) ~ p{E)E v l 2 . Hence flow prevents regions from 
getting stuck in progressively deeper traps and the aging 
process is truncated after a finite time. 

We are currently investigating more complicated non- 
linear strain histories [ p3| . In future work, explicit spatial 
structure and interactions between regions must be added 
so as to understand better the mutual dynamical evolu- 
tion of the attempt rate, the effective noise temperature 
and the disorder. One issue concerns the relative impor- 
tance of localized [|l9|,|8]-[33| versus avalanche-like [ p0|j2l| 
events in the relaxation of stress. 
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